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J3 ! Abstract 

We study, in the framework of open quantum systems, the geometric phase acquired by a uni- 

formly accelerated two-level atom undergoing nonunitary evolution due to its coupling to a bath 

of fluctuating vacuum electromagnetic fields in the multipolar scheme. We find that the phase 

variation due to the acceleration can be in principle observed via atomic interferometry between 
> 
ON ' the accelerated atom and the inertial one, thus providing an evidence of the Unruh effect. 
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When a quantum system undergoes a cyclic evolution, it may acquire a memory of this 
motion in the form of a geometric phase. This phase was first introduced by Pancharatnam 
while studying polarized beams passing through crystals [l]. In 1984, Berry studied the 
dynamics of a closed quantum system whose Hamiltonian varies adiabatically in a cyclic 
way, and found, besides the familiar dynamical phase, that there is an additional phase due 
to the geometry of the path enclosed during the evolution of the system in the parameter 
space [2]. Berry's work was soon generalized to nonadiabatic [3] and noncyclic evolution [4]. 
The geometric phase has so far been extensively studied, both theoretically and experimen- 
tally [5], and it has been fruitfully applied to many fields, such as the study of molecular 
dynamics [6] and electronic properties [7]. 

Recently, there has been interest in using the geometric phase for fault-tolerant quantum 
computation [8]. However, due to the inevitable interactions between the qubits and the 
environment, a pure state will be driven to a mixed state under the environment induced 
decoherence and dissipation. As a result, the geometric phase has to be generalized to general 
evolutions of an open system. Uhlmann was the first to define a mixed-state geometric phase 
via mathematical concept of purification [9]. Sjoqvist et al. put forward an alternative 
definition for the unitarily evolved nondegenerate mixed-state density matrix based on the 
interferometry [10]. This was soon generalized to degenerate mixed states by Singh et al. [11] 
and to the nonunitary evolution using the kinematic approach by Tong et al. [12]. Wang et 
al. defined a mixed-state geometric phase via mapping the density matrix to a nonunit vector 
ray in the complex projective Hilbert space [13]. Experiments based on NMR system [14] 
and single photon interferometry [15] have demonstrated the mixed-state geometric phase. 

As discussed above, the impact of environment on the geometric phase of open systems is 
an important issue in any practical implementations of quantum computing. In this regard, 
the effects of different kinds of decoherence sources on the geometric phase, such as dephasing 
and spontaneous decay, haven been analyzed [16]. In Ref [17], Rezakhani et al. have studied 
geometric phase for an open system, which is a spin-half particle in weak coupling to a 
thermal bath, and found that the phase varies with the temperature of the bath. Lombardo 
et al. [18] have studied not only how the geometric phase is modified by the presence of 



the different types of environments, but also estimated the corresponding times at which 
decoherence becomes effective. Chen et al. [19] focused on the geometric phase of an open 
two-level atom coupled to an environment with Lorentzian spectral density and explored the 
non-Markovian effect on the geometric phase. 

In quantum sense, every system, whatever it is, is an open system, since it is at least 
subjected to vacuum fluctuations. However, the geometric phase of an open system gener- 
ated by the nonunitary evolution due to its coupling to vacuum fluctuations is in general 
unobservable, as, practically, any phase variation is observed only via some kind of inter- 
ferometry between the involved state and certain selected reference states which are both 
inseparably coupled to vacuum. Nevertheless, if, somehow, vacuum fluctuations are modi- 
fied, then the geometric phase of the nonunitary evolution of an open system caused by its 
coupling to vacuum may become potentially observable. The modification of vacuum fluctu- 
ations induced by the acceleration of a two-level atom, for example, may provide such as a 
possibility, since, as is well-known, a uniformly accelerated observer perceives the Minkowski 
vacuum as a thermal bath of Rindler particles [20]. This is the so-called Unruh effect. So, 
the phase variation due to the acceleration of an two-level atom, which can in principle be 
observed through interference with an inertial atom, may provide evidence of the Unruh 
effect which is deeply related to the Hawking radiation. In this regard, let us note that many 
novel proposals have been suggested to detect the Unruh effect and the Hawking radiation 
in analog systems [21]. At this point, it may be worth pointing out that the Unruh effect is 
associated with quantization of the field in the Rindler accelerated frame. However, theoret- 
ical calculations performed from the perspectives of both the inertial frame and the Rindler 
accelerated frame with the Unruh thermal bath usually produce the same result on physical 
observables [22], as is the case in the weak decay of a uniformly accelerated proton [23], the 
bremsstrahlung effect associated with a uniformly accelerated point charge [24, 25], and the 
spontaneous excitation of a uniformly accelerated atom [26] 

Recently, Martin-Martinez et al. [27] have considered the possibility of using geometric 
phase to detect the Unruh effect. They examined an accelerated detector modeled by a 
harmonic oscillator which couples only to a single-mode of a scalar field in vacuum, and 



calculated the geometric phase acquired by the joint state of the detector and the field. As a 
result, cavities which are leaky to a finite number of modes are essential for the measurement 
of the acceleration influence in order to realize the single mode coupling and avoid the 
problems in the Unruh effect itself arising from introduction of boundaries. Such kind of 
cavities seems to be a major challenge in experimental implementation of their proposal. 
Here, we would like to consider a more realistic case and propose using the geometric phase 
of non-unitary evolution to detect the Unruh effect. We plan to study an accelerated two- 
level system which couples to all vacuum modes of electromagnetic (rather than scalar ) fields 
in a realistic multipolar coupling scheme [28]. We treat the accelerated two-level atom as an 
open system x in a reservoir of fluctuating vacuum electromagnetic fields and calculate the 
geometric phase of the accelerated open system undergoing non-unitary evolution because of 
the environment induced decoherence and dissipation. Since in our study, the atom couples 
to all vacuum modes, no cavity is needed in any experimental scheme to detect the phase. 
At this point, it is worth noting that the quantum geometric phase of an open system 
undergoing nonunitary evolution due to its coupling to a quantum critical bath has recently 
been demonstrated using a NMR quantum simulator [32]. 

Let us write the total Hamiltonian of the system (atom plus reservoir) as H = H s + 
H^ + H' . Here H s is the Hamiltonian of the atom, and, for simplicity, is taken to be 
H s — | Huqcts, in which <7 3 is the Pauli matrix, uq is the energy level spacing of the atom. H^ 
is the Hamiltonian of the free electromagnetic field, of which the details are not needed here. 
The Hamiltonian that describes the interaction between the atom and the electromagnetic 
field in the multipolar coupling scheme is given by H'{t) = —er ■ E(a;(r)) = —eJ2 mn r mn ■ 
E(x(r))cr mn , where e is the electron electric charge, er the atomic electric dipole moment, 
and E(x) the electric field strength. 

At the beginning, the whole system is characterized by the total density matrix p tot = 
p(0) ® |0)(0|, in which p(0) is the initial reduced density matrix of the atom, and |0) is the 



1 Let us note that the theory of open quantum system has been fruitfully applied to understand, from 
a different perspective, the Unruh, Hawking and Gibbons- Hawking effects, in Rcf. [29], [30] and [31], 
respectively. 



vacuum state of the field. In the frame of the atom, the evolution in the proper time r of 
the total density matrix p tot satisfies 

dptot(r) i 

dr =—jj[H,Ptot(T)]. (1) 

We assume that the interaction between the atom and the field is weak. In the limit of weak 
coupling, the evolution of the reduced density matrix p(r) can be written in the Kossakowski- 
Lindblad form [33, 34] 

^ = - Z j_[H eS ,p(T)]+£[p(T)}, (2) 

where 

1 3 

The matrix a^ and the effective Hamiltonian H e $ are determined by the Fourier and Hilbert 
transforms of the field correlation functions 

z 2 3 



G + (x -y) = : G -Y1 (+k*|-X-|ril+> (0\Ek{x)E j {y)\0) , (4) 



which are defined as follows 

g(\) = f" dre lXT G + (x(r)) , /C(A) = - f°° du ^\ . (5) 

J~oc ™ J_ 00 U- X 

Then the coefficients of the Kossakowski matrix a^ can be written as 

ay = A5ij - iBe ijk 5 k3 + C5 i3 5 j3 , (6) 

in which 

A = ±[g(uo) + G(-u>o)], B = ±[g(u>o)-g(-u>o)], C = -A. (7) 

The effective Hamiltonian H e g contains a correction term, the so-called Lamb shift, and one 
can show that it can be obtained by replacing uq in H s with a renormalized energy level 
spacing Q as follows 

1 h i 

H cS = -Mla 3 = -{u + -[/C(-w ) - £(w )]} °"3 • (8) 



For convenience, let us express the density matrix p in terms of the Pauli matrices, 
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(9) 



p(r) = Ul + J2 P i(r)aA 



Plugging Eq. (9) into Eq. (2) and assuming that the initial state of the atom is |^(0)) = 
cos ||+) + sin ||— ), we can easily work out the time- dependent reduced density matrix 

(-4At 2 fl I B-A ( -4At i\ 1 p -2(2A+C)r-iClT „:„ n \ 

2 + 2A l e 1J 2 e Smy (1Q) 

1 e -2(2A+C)r +4 nr sin Q l _ g -4Ar ^2 | _ B-A ^-4Ar _ ^ \ 

The geometric phase for a mixed state undergoing nonunitary evolution is given by [12] 

7 = arg NT v/A fc (O)A fc (T)(0 fc (O)|0 fc (T)) e "/o T ^MI^W)^j , (n) 

where Afe(r) and |0fc(r)) are the eigenvalues and eigenvectors of the reduced density matrix 
p(r). In order to get the geometric phase, we first calculate the eigenvalues of the density 
matrix (10) to get 

A ± (r) = ^(l±r / ), (12) 



where r) = yfp\ + e -^A+c)r sin 2 9 and p 3 = e~ AM co&B + f (e" 4Ar - 1). It is obvious that 
A_(0) = 0. As a result, the contribution comes only from the eigenvector corresponding to 

|0 + (r))=sin||+) + cos|e^|-), (13) 

where 



tan ^ = 7^. (14) 

The geometric phase can be calculated directly using Eq. (11) 

7 = -n / cos 2 — dr. (15) 

Jo 2 

Let us now calculate the geometric phase of an two-level atom which is uniformly ac- 
celerated, for example, in the x-direction. The trajectory of the atom is then described 

by 

2 
C (XT C CLT 

t(r) = - sinh — , x(r) = — cosh — , y(r) = z(r) = . (16) 

(X C (X c 



In order to get the explicit form of the geometric phase, we need the field correlation func- 
tions, which can be worked out using the two point function of the electric field 

(E M t))EMt>))) = ^(W - ^) | x _ xf _ (c 1 f _ cf ,_ g)2 • ("> 

The field correlation function for the trajectory (16) can then be evaluated from (17) in the 
frame of the atom to get 

G+(x x>) = e2 IH r l+>l 2 (m 

U [X ' X) 16^e hc 7 s i n h 4 if c (T-T'-i6)}- l J 

So, the Fourier transform of the field correlation function is 

„,^ A 3 e 2 |(-|r|+)|V a 2 \ ( , nc\\ , . 

°™ = Low l 1 + *» ) v + coth -J • (19) 

Consequently, the coefficients of the Kossakowski matrix ay and the effective level spacing 
of the atom are given by 



1 / a 2 \ e 27rca, °/ a + 1 1 
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*. = -ft = j7o ^1 + ^ ) ^^j , B.-^^ + Sj). (2°) 



where 70 = e 2 |(— |r|+)| 2 a;o/37r£o^c 3 is the spontaneous emission rate. Then the geometric 
phase can be obtained according to 

la = ~ - [ 1 ; ^a rfr , (22) 

7o 2 V Ve 4AoT sin 2 6 + (R - R e 4 ^ + cos fl) 2 / 

where i? = B a /A a . So, the phase accumulates as the system evolves, although the accumula- 
tion with time is not linear as in the unitary evolution case. For a single period of evolution, 
the result of this integral can be expressed as 

la = 5i [F(2n) - F(0)] , (23) 

where the function F(<p) is defined as 

1 1 . (l-Q 2 -R 2 + 2R 2 e 4A ^° _, . 
Ffr) = - 2 V-^{ 2R +S M 

--Ugn(Q) In (1 - Q 2 - R 2 + 2Q 2 e~ AA ^ + 2 \Q\3{<p) e~ 4A ^ M ) , (24) 



in which S((p) = ^R 2 e 8 ^^ + (1 - Q 2 - R^) e iA ^/^ + Q2 , Q = /? + C os# and sgn(Q) is 
the standard sign function. For small lo/oJo, which is generally true as we will see later, we 
can perform a series expansion to the result. For a single quasi-cycle, we find, to the first 
order 2 , 

7a « -tt(1 - cos 0) - it 2 — sin 2 0(1 + -%-= ] (2 + cos + -5 ^ cos ] . (25) 

The first term — tt(1 — cos 6) in the above equation is the geometric phase we would have 

obtained if the system were isolated from the environment, and the second term is the 

correction induced by the interaction between the accelerated atom and the environment. 

The geometric phase contains a term proportional to a 2 apart from the usual thermal term 

with a Planckian factor, and this term becomes appreciable when the acceleration is of the 

order of cuq, thus it enhances the accumulation of the geometric phase in contrast with the 

scalar field case where this term is absent. Let us note here that similar a 2 terms also appear 

in the studies of the energy shift [35] and the spontaneous excitation [36] of an accelerated 

atom once the scalar filed is replaced by the electromagnetic field. In the limit of a — > 0, 

which corresponds to the case of an inertial atom, there is still a correction, which comes 

from the zero point fluctuations of the Minkowski vacuum. The explicit form of this term 

reads 

7/M -vr(l -cos0) -tt 2 — (2 + cos 0) sin 2 . (26) 

2cu 

This correction is exactly the same as the one in Ref. [19], which is obtained by assuming an 
environment with a Lorentzian spectral density, and is very similar to the result in Ref. [37] 
derived from a different model. Thus the correction to the geometric phase purely due to 
the acceleration is 



$a = la ~ ll ~ -TT 



2 70 



a 2 ,_ „, / a 2 



(2 + cosfl)+ 1 + -— - - cosfl 



c 2 cu 2V V c 2 ul) e 2lTCU)a / a -I 



sin 2 9 . (27) 



2w 

This reveals that the geometric phase difference between the accelerated and inertial atoms 
depends on the properties of the atom (transition frequency Uq and the spontaneous emission 



2 Here we have omitted the Lamb shift terms, since it is obvious that these terms contain a factor 70/wo 
and they will only contribute to the phase at the second and higher orders of 70 /loq . 
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rate 70), the initial state (angle 9), and the acceleration a. If we assume that |(— |r|+)| is 
of the order of the Bohr radius ao, and uo of the order of E /h, where E Q = — e 2 /87r£ ao 
is the energy of the ground-state, then 70/^0 is of the order of 10~ 6 . For a given initial 
state, the phase difference increases with the acceleration, and it becomes significant when 
the acceleration is of the order cw . The initial state of the atom, i.e., the initial angle 9 in 
the Bloch sphere representation, also plays an important role. When 9 = and 9 = tt, which 
corresponds to an initial excited state and an ground state respectively, the phase difference 
vanishes, whereas it reaches its maximum in the regime near 9 = tt/2. For a typical transition 
frequency of the hydrogen atom, i.e., Uq ~ 10 15 s^ 1 , the acceleration needed to observe this 
effect is of the order of 10 23 m/s 2 , which is extremely high. However, if we consider two-level 
systems with lower frequency, the acceleration needed can be smaller. If we choose transition 
frequencies of the atom in the microwave regime, for example, uq = 2.0 x 10 9 s _1 , which is 
physically accessible [38, 39], then, for a = 4cu = 2.4 x 10 18 m/s 2 , the phase difference 
can reach 1.6 x 10~ 4 rad after a single period of evolution, which may be within the current 
experimental precision. 

The geometric phase discussed above, therefore, may be detected with an atom interfer- 
ometer. One first prepares the two-level atom in a superposition of upper and lower states 
in a Ramsey zone. In one arm of the interferometer the atoms move inertially, and in the 
other arm the atoms are accelerated. An interferometric measurement is taken when the 
atoms in the two arms meet. Here, let us recall that our calculations of the geometric 
phase are based on the comoving frame of reference of the atom. So, in the example we 
consider, according to Eq. (16), a single period for the accelerated atom in the comoving 
frame T = 2tt/uq ~ 3.1 x 10~ 9 s would transfer to a time interval of 5.1 s in the laboratory 
frame. Thus, one should prepare an inertial atom which moves fast enough so that a single 
period of time in its own frame also transfers to the same amount of time in the laboratory 
frame when the interference experiment is performed. A tricky point is whether the field to 
accelerate the atoms will change the structure of a real atom or even ionize it. 

Another delicate issue in experimental implementation is how to cancel the dynamical 
phase that the atoms acquire. For systems under nonunitary evolution like what we are con- 



sidering here, the removal of the dynamical phase from the total phase is a subtle issue [40]. 
However, since our purpose is to detect the Unruh effect associated with the acceleration 
of the atom, we do not really need a complete cancellation of the dynamical phase. In- 
stead, we may choose slightly different paths to control the relative dynamical phase to be 
much smaller than the geometric phase acquired in one period, so the result is effectively 
dominated by the geometric phase difference between the accelerated and inertial atoms. 
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